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Driven polymer translocation through a cylindrical nanochannel: Interplay between 

the channel length and the chain length 
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Using analytical techniques and Langevin dynamics simulations, we investigate the dynamics 
of polymer translocation through a nanochannel embedded in two dimensions under an applied 
external field. We examine the translocation time for various ratio of the channel length L to the 
polymer length A^. For short channels L -ti N, the translocation time r ~ N'^'^" under weak driving 
force F, while r ~ F^^L for long channels L ^ N, independent of the chain length A'^. Moreover, 
we observe a minimum of translocation time as a function of L/N for different driving forces and 
channel widths. These results are interpreted by the waiting time of a single segment. 

PACS numbers: 87.15.A-, 87.15.H- 



I. INTRODUCTION 

Polymer translocation through a nanopore or 
nanochannel is of considerable importance to a mul- 
titude of biological functions, such as protein and 
DNA transport through membrane channels and viral 
DNA injection into host cells [1|. Recent experiments 
show that monitoring single molecule translocation 
through protein or synthetic nanopores has potential 
applications in single DNA sequencing, single molecular 
characterization, and single molecular detection |2|-|23|. 

The transport of biopolymers through a nanopore has 
attracted broad interest in the statistical physics com- 
munity, as it rep resents a challenging problem in poly- 
mer physics [2J-l4J|. A quantity of particular interest 
is the average translocation time r as a function of the 
chain length N, usually assumed to follow a scaling law 
T ~ N". The scaling exponent a hereby reflects the 
efficiency of the translocation process. 

Standard equilibrium Kramers analysis of diffusion 
through an entropic barrier yields r ~ iV^ for the un- 
biased translocation and t ^ N for the driven translo- 
cation (assuming friction is independent of N) |24l. ua^. 
However, as Chuang et al. [23] noted, the quadratic scal- 
ing behavior for unbiased translocation cannot be cor- 
rect for a self-avoiding polymer. The reason is that the 
translocation time is shorter than the Rouse equilibra- 
tion time of a self-avoiding polymer, Tfi ^ N^^'^'^, where 
the Flory exponent ly = 0.588 in three dimensions and 
ly = 0.75 in two dimensions |45|, thus rendering the con- 
cept of equilibrium entropy and the ensuing entropic bar- 
rier inappropriate for translocation dynamics. Chuangei 
al. [27[ performed numerical simulations with Rouse dy- 
namics for a two dimensional lattice model to study the 
translocation for both phantom and self-avoiding poly- 
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mers. Their results show that for large N translocation 
time T ~ A^^+^'^j which scales approximately in the same 
manner as the equilibration time but with a much larger 
pref actor. For driven translocation, Kantor and Kardar 
[3l| have demonstrated that the assumption of equilib- 
rium in polymer dynamics by Sung and Park |2J] and 
Muthukumar [25[ breaks down more easily and provided 
a lower bound r ~ N^'^'^ for the translocation time by 
comparison to the unimpeded motion of the polymer. 
Most recently, we have found that for faster transloca- 
tion processes a — 1.37 in three dimensions [46[, while it 
crosses over to a = 1 -|- i' for slower translocation, cor- 
responding to weak driving forces and/or high friction 

M- 

Translocation process can be affected by many factors, 
such as the driving force, and the shape and the size of 
the pore. As to the driving force, it can be provided by 
an electric field, a chemical potential, a pulling force, or 
geometrical confinement of the polymer. 

However, most of previous studies focus on short 
nanopores, where the pore length is much smaller than 
the radius of gyration of the polymer. The dynamics 
of polymer translocation through a long nanochannel is 
of great importance, because it is related to many tech- 
nological applications, such as the ultrafiltration process 
[47| . Although polymer translocation through a cylindri- 
cal channel of finite diameter and length between two 
spherical compartments has been theoretically investi- 
gated [48|, |49| , equilibrium entropy of polymer as a func- 
tion of the position during the translocation is used, lead- 
ing to the results to be correct at most for Gaussian chain. 
The translocation dynamics for polymer through a long 
channel is still not clear due to the interplay of the chain 
length and the length of the channel. To this end, us- 
ing both analytical techniques and Langevin dynamics 
(LD) simulations, we investigate the dynamics of poly- 
mer translocation through a long nanochannel. 

The paper is organized as follows: in Section H, we de- 
scribe our simulation model and methods; in Section HI, 
we present the results of blob and scaling theories with 




conservative, frictional, and random forces, respectively, 
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FIG. 1: Schematic representation of polymer translocation 
through a two-dimensional long nanochannel of length La 
and width Da. The driving force acting on the bead in the 
channel can be provided by an applied intra-channel electric 
field. For long nanochannels, the translocation process can 
be broken into three components r ~ ti + T2 -(- T3, where 
Ti, T2 and rs correspond to initial filling of the nanochannel, 
transfer of the polymer from the starting cis side to the trans 
side, and finally emptying of the nanochannel, respectively. 



discussing our simulation data of polymer chain translo- 
cation process; the conclusions are drawn in Section IV. 



II. MODEL AND METHODS 

We consider the two dimensional geometry where the 
nanochannel is sandwiched between two parallel walls 
that are close to each other with the polymer chain, 
as shown in Fig. [TJ In our simulations, the polymer 
chains are modeled as bead-spring chains of Lennard- 
Jones (LJ) particles with the Finite Extension Nonlinear 
Elastic (FENE) potential [50|. Excluded volume inter- 
actions between beads are taken into account by a short 
range repulsive LJ potential 



ULjir) 



(a/r)^^ - [a/rf 



r > 2I/V, 



(1) 



where a is the diameter of a bead and e is the potential 
depth. The connectivity between neighboring beads is 
modeled as a FENE spring with 



UpENEir) 
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fci?2 In (1 - rVi?2) 



(2) 



where r is the distance between consecutive beads, k is 
the spring constant, and Rq is the maximum allowed sep- 
aration between connected beads. 

As shown in Fig. [TJ the wall is formed by stationary 
particles with the diameter a. In unit of cr, the nanochan- 
nel has length L and width D. Between all bead-wall 
particle pairs, there exists the same short range repul- 
sive LJ interaction as described above. In the Langcvin 
dynamics (LD) simulations, each bead is subjected to 



mfj = -V([/lj + Ufene) 



Cv, 



(3) 



where m is the bead mass, ^ is the friction coefficient for 
a single bead, v^ = f j is the bead velocity, and Ff is the 
random force satisfying the fluctuation-dissipation theo- 
rem. The external force is expressed as Fcxt = Fx, where 
F is the strength of force exerted on the translocating 
beads located inside the channel, and i is a unit vector 
in the direction along the channel. Experimentally, this 
external driving force acting on the bead in the channel 
can be provided by an applied transmembrane electric 
field. 

In the present work, we use the LJ parameters s and 
a and the bead mass m to fix the energy, length and 
mass scales, respectively. The time scale is then given 
by ^Lj = (nicr^/e)^'^. The dimcnsionlcss parameters in 
our simulations are ^ = 0.7 and fc^T =1.2 unless oth- 
erwise stated. For dimensionless parameters for FENE 
potential, we use Rq = 2 and k = 7 [54|. We have also 
checked that using Rq = 1.5 and fc = 30 does not change 
our results. The Langevin equation is integrated in time 
by a method described by Ermak and Buckholz [5J] in 
two dimensions. 

Initially, the first bead of the chain is placed just in- 
side the channel (at x = 0.75, y = 0), while the remain- 
ing beads are under thermal collisions described by the 
Langevin thermostat to obtain an equilibrium configu- 
ration. The translocation time is defined as the time 
interval between the entrance of the first segment into 
the pore and the exit of the last segment. Typically, we 
average the data over 1000 independent runs. 



III. RESULT AND DISCUSSION 

A. Theory 

For translocation through short channels, the translo- 
cation time T has been estimated by Kantor et al. |3l| 
using the scaling arguments, where an essential assump- 
tion is that the chain is not severely deformed during 
translocation and in particular the chain configurations 
at both cis and trans sides are close to equilibrium. For 
slow translocation, these assumptions are satisfied. Thus, 
the translocation time can be written as |3l| 

T - Rg/v - N^+'^/F (4) 

for L <^ Rg. Here, L is the channel length, Rg ^ N" is 
the radius of gyration of the polymer without the con- 
finement, and V ~ F/{N£^) is the average translocation 
velocity. The predicted exponent 1 + 1/ is observed in the 
slow dynamics regime [46[ . 

For translocation through long channels, we can break 
down the translocation process into three components, 



10^ 



10' 



F=0.25 



Long channel regime 
Slope= 1.0 



Short channel regime 
Slope= -0.75 




0.01 




L/N 



100 



FIG. 2: t/N as a function of L/N for the channel width Z) = 3 
and the driving force F = 0.25. Here, we use A*' = 32, 64, 96 
and 128. 
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FIG. 3: Waiting time distribution for TV = 50, D = 3, F = 
0.25 and different cliannel length L. Here, ta=i includes the 
the time duration for initial filling process ri. 



as shown in Fig. [T] The total translocation time can be 
written as a sum of three contributions r ss ti + T2 + T3, 
where ri, T2 and T3 correspond to initial filling of the 
nanochannel, transfer of the polymer from the starting 
CIS side to the trans side, and final emptying of the 
nanochannel, respectively. Particularly, for an infinite 
long channel of length L, the translocation dynamics is 
dominated by r2, which can be written as T2 ^ L/v, 
with V being the average translocation velocity. Based on 
the balance of the frictional force and the driving force, 
NS,v — NF where F is the driving force acting on one 
bead, we have v — F/^. Therefore, the translocation 
time 



r « T2 - L/f - L^/F 



(5) 



for L ^ Rg. This result indicates that t is independent 
of the chain length N. 

Based on Eqs. ^ and (O, there exists a crossover 
for the scaling behavior of r as a function of L. For 
short channels L <^ N, the translocation time r ^ N^^'^ 
under weak driving force F [4y], while r ^ F^^L for long 
channels L ^ L. However, r is still not clear for various 
scaled pore length L/N. To this end, in the following we 
give more details from simulations. 



B. Numerical results 

In our simulations, we find that the translocation prob- 
ability is approximately independent of the chain length 
TV as observed in the previous study ^^. This result is 
also in agreement with the experiments |54l . |55| , where 
the polymer chain insertion process is entirely local to 
the portion of the chain at the interface and independent 
of the chain length. 

A particularly interesting question concerns r for var- 
ious L/N. Fig. [5] shows the scaled translocation time 





□ -□-(ts=i+ts=2+-+ts=N)/N 


10- 


\ -^-(ts=2^ts=3^-*ts=N)/N 


8- 
0) 




E 6- 

4-1 


\ uy^^=^ 


4- 


^^^>^^£^2l^ 


2- 


y^'"'^^"^-''''-^ 


0- 





0.0 



0.2 



0.4 



0.6 



0.8 

L/N 



1.0 



1.2 



FIG. 4: Sum of waiting times for all beads {^^^^ts^i/N), 
sum of waiting times for beads numbered from 2 to TV 
(5^;^2 *s=i/A'') and waiting time for the first bead ts=i/N 
as a function of L/N for iV = 50, D = 3 and F = 0.25. Here, 
is=i includes the the time duration for initial filling process 

Tl. 



t/N as a function of L/N for F = 0.25. For L/N < 0.1, 
we find t/N ^ {L/N)-" with v = 0.75 being the Flory 
exponent in two dimensions. This implies t '^ N°' with 
a = 1 -|- I/, in agreement with the prediction in Eq. ^ 
for L<^Rg. 

In our previous results, we find that with increasing the 
driving force F the exponent a crosses over from 1 -I- 1/ to 
2j/ in two dimensions [32| and from 1 + t/ to 1.37 in three 
dimensions |46| . Due to translocation through a quite 
short pore with a small L, the driving force acting on 
the chain is quite weak, and thus the exponent 1 + i^ is 
in agreement with the previous result for slow dynamics 
regime. However, in contrast to our results |32l . 14211 . re- 
cently Lehtola et al. [501 and Dubbeldam et al. [531 ^^^o 
observe the crossover exponent from 2z^ to 1 -\- v with 



increasing F. This discrepancy may be from the chain 
length and the driving force used in the simulations |58| . 

When L/N > 10, we observed t/N - L/N, which 
demonstrates that t ^ L, independent of the chain 
length. This scaling behavior is in agreement with the 
prediction in Eq. (O for L ;» i?g. Therefore, the ob- 
served crossover scaling behavior is in good agreement 
with our theoretical predictions in Eqs. @ and ([5]). 

The most interesting result in Fig. [2]is that the translo- 
cation time has a minimum for an optimal value of L/N 
at L/N « 0.38. To understand the existence of the 
minimum, we take into account the dynamics of a sin- 
gle segment passing through the channel during translo- 
cation. We numerically calculated the waiting times 
for all beads in a chain. The waiting time of bead s 
for successful translocation hereby is defined as the av- 
erage time between the events that bead s and bead 
s + 1 exit the channel. Fig. [3] shows the waiting time 
distribution for iV = 50, Z? = 3, F = 0.25 and dif- 
ferent channel length L. For L > 5, it takes much 
longer time for the first segment to exit the channel, 
denoted as ts=i- Here, we should point out that ts=i 
includes the time duration for initial filling process ti 
because the first segment must pass through the whole 
channel before its exiting the channel. Moreover, ts^i 
increases with increasing L, while the waiting times for 
other monomers (is=2, is=3---^s=Ar) initially decrease and 
then saturate with increasing L. Quantitatively, we mea- 
sure the sum of waiting times for all beads J2i=i ts=i/N, 
sum of waiting times for beads numbered from 2 to N, 
'Ylii=2^s=i/N , and the waiting time for the first bead 
ts=i/N as a function of L/N, see Fig. HI With increasing 
L/N, ts=i/N first increases rapidly followed by a slower 
increases, while X]i=2 ts=i/N initially decreases rapidly 
and then almost saturates for long channels. The inter- 
play between ts=i/N and X]i=2 ts=i/N leads to the ob- 
served minimum of the translocation time (the curve for 
Y^i=i'^s=i/N). Recently, this effect was also observed in 
unbiased translocation through long channel based on a 
computer simulation [28|, and was predicted for translo- 
cation through a long channel between two spherical com- 
partments [49[. 

To generalize above observed results, we examine t/N 
as a function of L/N for different driving forces and chan- 
nel widths in the following. 



L The influence of the driving force on the translocation 
time 

Fig. [5] shows t/N as a function of L/N for N = 50, 
D — i and different driving forces. As expected, with 
increasing F, t/N decreases for the range of L/N we 
examined. Moreover, we observe the minimum of the 
translocation time, {t /N)mim for different F. The crit- 
ical ratio, {L/Na)c, corresponding to the minimum of 
t/N , almost does not change with F. 
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FIG. 5: t/N as a function of L/N for TV = 50, D = 3 and 
different driving forces. 
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FIG. 6: The minimum of the translocation time as a function 
of the driving force for A*' = 50 and _D = 3. 
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FIG. 7: Total translocation time r as a function of the driving 
force F for different L/N. The chain length N = 32. 




FIG. 8: T /N as a function of L/N under the driving force 
F = 0.25 for the chain length Af = 50 and different channel 
widths D. 



As shown in Fig. [51 wc further find {T/N)min ^ F^ ^ 
with /3 — —0.90 ±0.01. This exponent is close to -1. 
Fig. [7]shows the translocation time r as a function of the 
driving force F for different L/N . r initially decreases 
and then increases with increasing L/N ^ which further 
confirms that r has minima at L/N k, 0.4 for different 
driving forces F. In addition, we also observe r ^ F~^ 
for long channels with L/N > 2.0. For L/N < 0.4, 
exponents are a little larger than -1. Particularly, for 
L/N = 0.0625 and F < 0.2, the driving force is very 
weak and the translocation is almost controlled by pure 
diffusion as demonstrated in Ref. |4J, leading to a weaker 
dependence of r on F. 



The influence of the channel width on the translocation 
time 




FIG. 9: The minimum of the translocation time as a function 
of the channel width D under the driving force F — 0.25 for 
AT = 50. 
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FIG. 10: t/N as a function of the channel width D ior N = 50 
and F = 0.25. 



t/N as a function of L/N under the driving force 
F = 0.25 for iV = 50 and different D is shown in Fig. 
|S1 where r also shows the minimum. The critical ra- 
tio, {L/N)c, corresponding to the minimum of r/7V, only 
has a very weak D dependence based on the fitted curves. 
This behavior is a little different from the results for non- 
driven translocation through a long channel between two 
spherical compartments j49| . where the critical channel 
length corresponding to the minimum of r always in- 
creases with increasing the channel width. This is be- 
cause the entropic force for polymer partially inside the 
channel depends on D for unbiased case. 

In addition, t/N decreases with increasing D for dif- 
ferent L/N. This is because the channel can accom- 
modate more beads with increasing D, leading to larger 
driving force. The minimum of the translocation time, 
{t /N)min decrease rapidly first, and then is slowed down 
with increasing D, see Fig. [HI For _D > 3, we ob- 
serve {T/N)rmn ^ D-O^OiO-Ol _ ^Ig. [TU] showS thc t/N 



as a function of thc channel width D for A^ = 50 and 
F ~ 0.25. t/N has a weaker dependence on D compared 
with the equilibrium longitudinal size of the chain con- 
fined to a channel, R\\ - Na{%Y/''-^-^ ~ ND'^/^ H. 
This may be due to the elongations of the chain inside the 
channel under the driving force and the non-equilibrium 
process of the translocation. 



IV. CONCLUSIONS 

Using analytical techniques and Langevin dynamics 
simulations, we investigate the dynamics of polymer 
translocation through a nanochannel embedded in two 
dimensions under an applied external field. We examine 
the translocation time for various ratio of the channel 
length L to the contour length of polymer TV. For short 
channels L <^ N , the translocation time r ^ iV^"*"" under 



weak driving force F, while r '^ F~^L for long channels 
L ^ N ^ independent of the chain length N . Moreover, 
we observe a minimum of translocation time as a function 
of L/N for different driving forces and channel widths. 
These results are interpreted by the waiting time of a 
single segment. This work is relevant for the applied 
biophysics field. Our results can help design a device ge- 
ometry to achieve higher separation resolution for DNA 
molecules using nanochannel. 

In our previous study [4j], we have examined the dif- 
ferent case where the driving force acts solely on the bead 
at the entrance of the channel. Due to the crowding ef- 
fect induced by the partially translocated monomers, the 



translocation dynamics is significantly altered in compar- 
ison to an unconfined environment. The scaling exponent 
of the translocation time r with the chain length N de- 
pends on the channel width and the driving force. 
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